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ABSTRACT 

An investigation of helicity injection by photospheric shear motions is presented for two active 
regions, NOAA 11158 and 11166, using line-of-sight magnetic field observations obtained from 
the Helioseismic and magnetic Imager on-board Solar Dynamics Observatory. We derived the 
horizontal velocities in the active regions from Differential Affine Velocity Estimator (DAVE) 
technique. During the six day evolution period of the active regions, we found persistent strong 
shear motions at the maximum velocity in the range of 0.5-0.7 km s^along the magnetic polarity 
inversion line and outward flows from the peripheral regions of the sunspots. The helicities 
injected in NOAA 11158 and 11166 during the six days' period were estimated as 13.30 x 10 42 Mx 2 
and 9.5x 10 42 Mx 2 , respectively. Temporal profiles of helicity injection showed impulsive variations 
at the onset times of flares/CMEs due to the negative helicity injection in the dominant region 
of positive helicity density. The spatial examination of helicity density maps showed that these 
variations resulted mainly due to the negative helicity injection in the regions of opposite helicity 
that were co-spatial with flaring sites. These co-spatial and co-temporal variations of helicity 
injection with flares are interpreted to be due to the relaxation from the state of high shear by 
the observed motions to a lower or s hear-free state r eleasi ng energy in the form of eruptive events. 
This agrees with the simulations bv lKusano et al. (|2004 ) for the triggering mechanism of flares. 
However, for the flares of smaller magnitude no clear evidence of such changes was available. We 
expect that this difficulty could be addressed by careful calculations for the emergence term of 
helicity injection using the vector magnetic field data once they are available. Our study suggests 
that the existence of opposite helicity fluxes can trigger eruptive events, viz., flares and CMEs, 
and promises to be useful in forecasting the transient activity of ARs. 

Subject headings: Sun: activity — Sun: flares — Sun: magnetic fields — Sun:Coronal Mass ejections 



1. Introduction 

Magnetic helicity is an important topological 
property of solar active regions (ARs) and is a 
measure o f its content, viz., twist and writhe of the 
field lines (jBerger fe Fieldlll984trFinn fc Antonsenl 
19851 ) . It is gauge invariant for a closed volume of 
space. Sun's outer atmosphere is dominated by 
magnetic field at all scales. Dynamic phenomena, 
such as, energetic flares and coronal mass ejections 
(CMEs) occur due to the loss of equilibrium dur- 



ing the evolution of magnetic fields in solar ARs. 
Magnetic helicity has become an important phys- 
ical parameter in the context of solar transient 
phenomena. It is one of the few global quanti- 
ties which is conserved even in resistive magneto- 
hydrodynamics on a timescale less than that of the 
global diffusion. There exists no absolute measure 
of helicity within a sub- volume of space if that 
sub- volume is not bounded by a magnetic surface. 
However, a topologically meaningful and gauge in- 
variant relative helicity for such volumes can be 
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measured. 

There are several methods for estimating helic- 
ity in solar ARs. By the force- free field assumption 
of coronal magnetic field, we have: 



V B = aB 



(1) 



where a is the force-free parameter, also known 
as helicity or twist parameter. Assuming a to 
be constant for the whole AR, we can fit ob- 
ser yed vector magnetograms to deduce the value of 
a (IPevtsov et all 19951: lHagvard fe Pevtsovl fl999: 



Tiwari et al.l I2009T ) . Latitudinal variation of he- 



licity of photospheric magnetic fields, statistical 
significance of the observed temporal variations of 
the ARs' hemispheric helicity rule, as measured 
by the latitudinal gradient of the best-fit linear 
for ce-free parameter a. etc., have been discussed 
bv lPevtsov et al.l ((20081 ) . 

The Pointing theore m for helicity in an op en 
volume as derived by iBerger fe Field! ( 1984 ) is 
given by: 



= j 2(B t ■ A p )v z ds - j> 2(A p • w)B z ds (2) 

where A p is the vector potential of the magnetic 
field, B, which is uniquely specified by the ob- 
served flux distribution on the surface as 



V x A p • h 



B x 



V-A p = 0; A E 



(3) 



Equation [5] shows that the helicity of magnetic 
fields in an open volume may change by the pas- 
sage of helical field lines through the surface (first 
term) and/or by photospheric foot point motions 
of the field lines (second term). The temporal 
evolution of magnetic helicity on the photosphere 
characterizes both the injection of magnetic he- 
licity from the convection zone into the solar at- 
mosphere and the changes in the coronal magnetic 
field configurations related to eruptive events, such 
as the CMEs, propagating into the interplanetary 
medium. 

During the past years, several attempts have 
been made to estimate magnet i c helic ity from suit- 
able solar observations. IChael (|200ll ) developed a 
method for determining the helicity flux (the sec- 
ond term in Equation [5]) passing through the pho- 
tosphere. They used a time series of photospheric 



line-of-sight (LOS) magnetograms to determine 
horizontal velociti es by local correlation tra cking 
(LCT) technique ([November fe Simonl Il988t ) . Us- 
ing this method, vector potential A p was con- 
structed by using photospheric B z as boundary 
conditions with Coulomb gauge in terms of Fourier 
Transform (FT) as: 



A px = FT- 1 



A pv = FT- 1 



jky 



L.2 _i_ L,2 

Z?2 _i_ 1,2 
K x + K y 



FT (B z ) 



FT (B z ) 



where k x , k y are spatial frequencies in the x, 
y directions, respectively. Later, this method 
was applied to many ARs by several authors 
(Chae et alj|200lt iMoon et alj|200 j; iNindos et al' 



2003 [IChae et al.l2004HWelsch et al.|2004UJeong fc Chae 
20071 ). However. IPariat et al.l ([20051) showed that 
this method of calculation introduced artificial po- 
larities of both signs in the helicity density maps 
with most of the flow patterns. Hence, they sug- 
gested to use relative velocities for calculating the 
helicity injection rate: 



dH 
~dt 



-1 
2^ 



[(x 



S JS' 



X(x')5„(x)dS'dS 

(4) 



where u is the foot-point velocity at the position 
vector x, and B n is the normal component of the 
observed magnetic field. This equation shows that 
the helicity injection rate can be understood as 
the summation of rotation rates of all the pairs 
of elementary fluxes weighted with their magnetic 

flux. 

Furthermore, ISchuckl ([2005h have shown that 
the LCT method is inconsistent with the magnetic 
induction equation, which governs the tempo- 
ral evolution of the photospheric magnetic fields. 
Tracking methods have serious limitations as they 
largely underestimate the a mount of helicity in- 
jected by the shear term (|Demoulin fc Bergerl 
2003). They reported that the magnetic energy 
and helicity fluxes should be computed only from 
the horizontal motions deduced by tracking the 
photospheric cross-section of magnetic flux tubes. 
These horizontal motions include the effect of both 
the emergence and the shearing motions. They 
analyzed the observational difficulties involved in 



2 



deriving such fluxes, in particular, the limitations 
of the correlation tracking methods. One of the 
main limitations has been the coarse spatial res- 
olution of the available observations which limits 
the deduced velocities to the velocity correspond- 
ing to the group motion of unresolved bunch of 
thin flux tubes covered by a pixel. Also, tracking- 
motions have difficulties in the areas where suffi- 
cient contrast is lacking, such as in the sunspot 
umbrae. 

Several alternative, improved methods have 
been developed for inferring plasma velocities 
consistent with the induction equation at the 
photospheric level, based on the LOS, as well 
as, vector magnetogram s. Induction method 
(IM lKusano et al.1 d2002hh induction local corre- 
lation trackingQLCT rWelsch etaH fa004) . min- 
imum energy fit (MEF; iLongcopd (|2004)), dif- 
ferent i al affi ne velocity estimator (DAVE; SchuclJ 
(120051. 12006^ and differential affine velocity es- 
timator for vector magnetograms (DAVE4VM; 
Schuckl (120081) ) have been developed for the deter- 
mination of horizontal component of motion. On 
the other hand, the vertical component of velocity 
can be determined from the doppler observations 
near the disk center. DAVE4VM method on the 
other hand requires vector magnetograms. The 
performa nce of different tech niques have been ex- 
amined in lWelsch et al. ( 2007 ') which showed that 
the MEF, DAVE, FLCT, IM, and ILCT algo- 
rithms performed comparably. Furthermore, they 
reported that while the DAVE estimated the mag- 
nitude and direction of velocities slightly more ac- 
curately than the other methods, MEF's estimates 
of the fluxes of magnetic energy and helicity were 
more accurate. 

Time series data of photospheric magnetograms 
have been extensively used to derive magnetic 
helicity and its evolution in order to examine 
its ro le in the level of transient activity of the 
ARs. |m 

oon et al.l ([20021 ) reported impulsive varia- 
tions of hel i city dur ing some M and X-class flares. 
Park et al. I (|2010bl) have also studied the solar 
flare productivity in relation to the helicity in- 
jection using a large sample of 378 active region. 
Using SOHO-MDI magnetograms, they reported 
variation of helicity injection rates and a signifi- 
cant helicity accumulation of (3 — 45 x 10 42 ) Mx 2 
over several days around the time of flares above 
M5.0. Most of the previous studies that used data 



from Michelson Doppler Imager (MDI) onboard 
SOHO had the time resolution of 96 minutes. The 
rather coarse time resolution between two consec- 
utive observations has been a matter of concern 
in the above calculations because the contribution 
from the motion of short lived magnetic features 
in small intervals is difficult to be accounted suit- 
ably. This underlines the need for observations of 
magnetic fields with high temporal resolution. 

The above mentioned issues can now be ad- 
dressed with the availability of a better cadence 
of 12 minutes by the recently launched Helioseis- 
mic and magnetic Imager (HMI) onboard Solar 
Dynamics Observatory (SDO). Our main objec- 
tive in the present work is to reinvestigate the 
role of helicity injection in relation with flares 
and CMEs using the high-resolution data obtained 
from SDO-HMI. We intend to utilize this oppor- 
tunity to revisit some of the previous studies in- 
volving computations of helicity rate for two ARs, 
NOAA 11158 and NOAA 11166, that appeared 
during February and March 2011, respectively, in 
the ascending phase of the current Solar Cycle 24. 

Using the high quality HMI data obtained for 
the two ARs, we intend to examine whethe r the 
variations as reported by iMoon et al. (2002) and 
Park et al. (2010) for energetic flares also oc- 
curred during the flares of lower magnitude. It is 
of particular interest to investigate if such changes 
were associated with the CMEs as well. We have 
used DAVE technique for retrieving horizontal 
foot-point velocity from the LOS magnetograms. 
Thereafter, using Equation |4] we have determined 
helicity injection rates and the accumulated helic- 
ity in the two ARs due to foot-point shearing mo- 
tions during their disk transit. It has been inferred 
from the previous studies that most of the helicity 
injection corresponds to magnetic flux emergence 
in the ARs. We, therefore, attempt to interpret 
the variations found in these physical parameters 
in relation to the occurrence of flares and CMEs. 
In particular, we investigate whether the transport 
of magnetic helicity plays a role in solar eruptions. 

We organize this paper as follows. Description 
of the data used in this study and the procedures 
of the data processing are given in Section [2] Re- 
sults obtained for the two selected ARs are pre- 
sented in Section [3] and the discussions are given 
in Section [4] The summary and conclusions are 
given in Section [5j 
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2. Data and Methods of Analysis 

For our study, we have used high resolution 
LOS magnetograms at a cadence of 12 minutes 
obtaine d from the Heliosei smic Magnetic Imager 
(HMI; ISchou et al.1 (1201 lb ) onboard Solar Dy- 
namic Observatory (SDO). HMI observes the full 
solar disk in the Fe I 6 173 A spectral line with a 
spatial resolution of 1 arc-second. HMI provides 
four main types of data: dopplergrams (maps of 
solar surface LOS velocity), continuum filtergrams 
(broad- wavelength maps of the solar photosphere), 
LOS and vector magnetograms (maps of the pho- 
tospheric magnetic field). 

NOAA 11158 (19°S) and 11166 (10°N) ap- 
peared on the disk during February 11-20, 2011 
and March 03-16 2011 respectively. These ARs 
were very active, and produced some intense X- 
class flares associated with CMEs in addition to 
many M- and C- class flares during their disk 
transits. From the AIA observations, intermittent 
mass expulsions were seen, many of which turned 
into large, fast moving CMEs further confirmed 
by STEREO. Table [1] gives a list of flares^ (as 
recorded by GOES) and CMEfB 

Magnetograms obtained at different times were 
aligned by the method of Chae et al. (2004). In 



this method, an image of the AR taken at the cen- 
tral meridian is considered as the reference image. 
All other images, in time acco unted for differential 
rotation ( Howard et al. 1990f) along with the lati- 
tudinal difference of center of reference image from 
the ephemeris information, were remapped on to 
the disk center. This method is efficient in reduc- 
ing errors due to geometrical foreshortening and 
the AR is transformed to the disk center. The dif- 
ference between the normal and LOS component 
was corrected by cosine of the distance of the AR 
center from the disk center by assuming the hori- 
zontal field co ntribution for the transformatio n to 
be negligible ( Venkatakrishnan fe Gary 19891 ). 

We followed the transits of the two selected ARs 
from east to west on the solar disk. In order to 
have negligible errors in geometric correction, we 
restricted ourselves to a region within ±40°. With 



1 Obtained from the website www.solarmonitor.org 

2 by scrutinizing AIA 304 A quicklook images mirrored at 
http: //jsoc . st anf ord. edu/dat a/ ai a/ images/2011/ 

and further confirmed by the timings from 
http : //spaceweather . gmu . edu/seeds/ 



this constraint, we confined our study of the tem- 
poral evolution of the ARs to six days's period 
around their central meridian passage. 

We derived the horizontal velocities due to 
shearing motion of foot-poi nts on the pho tosphere 
by using DAVE technique ISchuckl ([20061 ) ). The 
DAVE technique is essentially a local optical flow 
method that determines the mass velocities within 
the windowed region. Further, it adopts an affinc 
velocity profile specifying velocity field in the win- 
dowed region about a point and constrai ns that 
profile to satisfy the induction equation. ISchuck 



(|2008l ) found the best performance of this method 
at approximately a square window of pixels. Any 
tracking method depends on two parameters i.e, 
the window size and the time interval. For a 
given time interval At, the window should be large 
enough so that tracked features remain within the 
window. Also, it should be small enough to vali- 
date locally with affine velocity profile. Since the 
ARs were evolving rapidly, we chose a window size 
of 21x18 pixels after a careful verification of the 
physical flux motions. The dependance of helicity 
injection rate on window size and time difference 
between the tracked maps using this method were 
verified. We have taken velocities greater than 0.7 
km s _1 as unphysical, and hence, set to zero in 
our study. Moreover, as the HMI magnetic field 
measurement precision is 10G, we have set this as 
the threshold to avoid errors while retrieving ve- 
locities. Further d etails of thi s meth od are given 

(l201lh . 



in a recent work of Tian et al 



Computation of t he helicity r ate u sing the 
method proposed by iPariat et al.l (|2005 ) at each 
pixel of the AR map (c.f., Equation!!]) is a tedious, 
time consuming process. However, we chose to use 
this method for reducing the effect of fake polar- 
ities of helicity. By restricting the calculations 
at pixels above the threshold (>10G) helps in re- 
ducing the computation time typically by 15-25%. 
Parallelization in integrand computation further 
reduces the time approximately by a factor of the 
number of proc essors used. The same equation 
as rewritten by IChaei ( 2007 ) to suit the convolu- 
tion algorithm by Fourier transform is faster than 
the direct integration method. The intrinsic prob- 
lem of Fourier transform with periodicity could 
be overcome by padding the array correspond- 
ing to the data points to get results as obtained 
by direct integration method. In this study, we 
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Table 1: List of Flares and CMEs 



AR 


Date 


Flares 


CMEs 


(NOAA) 


dd/mm/yyyy 


magnitude (time) 


(time) 


11158 


11/02/2011 


No flares 


No CMEs 




12/02/2011 


No flares 


No CMEs 




13/02/2011 


C4.7(13:44),C1.1(21:17),M6.6(17:28) 


21:30,23:30 




14/02/2011 


C8.3(04:29),C7.1(03:22), CI. 8(08:39) 
C7.0(13:47),C1.2(23:14),C2. 7(23:40) 


02:40,07:00,12:50,17:30,19:20 
14:00,15:00,17:30,18:40,19:30 




15/02/2011 


C1.0(10:02),X2.2(01:44),C2. 7(00:31) 
C4.8(04:27),C4.8(14:32), CI. 7(18:07) 
C6.6(19:30),C1.3(22:49) 


00:40,02:00,03:00,04:30,05:00 
09:00 




16/02/2011 


C2.0(00:58),C5.9(05:40),C2. 2(06:18) 
C9.9(09:02),C3.2(10:25),C1.0(11:58) 
M1.0(01:32),M1.1(07:35),M1. 6(14:19) 


14:35 


11166 


06/03/2011 


C3.9(15:21) 


02:00,15:20 




07/03/2011 


Ml. 7(13:45) 


14:25,22:10 




08/03/2011 


C7.7(23:10) 


14:30,19:00 




09/03/2011 


C7.3(10:35),X1.5(23:13) 


06:40,10:40,21:45 




10/03/2011 


C2.9(03:50),C4.7(13:42),C2.0(14:21) 
C4.0(19:00),M1. 1(22:34) 


04:50,07:10 




11/03/2011 


C1.4(00:29),C1. 1(01:46), C2. 8(02:24) 
C5.5(04:15),C4.3(07:22), CI. 1(08:13) 
C2.0(11:10),C3.6(11:43), CI. 1(16:04) 


00:50 



All times are in UT 
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have implemented the former to get sufficiently 
accurate results. 

3. Evolution of Magnetic Flux and Helic- 

ity 

The evolution of observed magnetic flux and 
the computed helicity rates are presented in the 
following for the two selected ARs NOAA 11158 
and NOAA 11166 with the methods and proce- 
dures explained in Section [2J 

3.1. AR NOAA 11158 

This AR appeared as small pores at the helio- 
graphic location E33S19 on 2011 February 11 as 
seen in the full disk HMI photohcliograms. There- 
after, it grew very rapidly during the next two 
days as the small pores merged and formed big- 
ger sunspots. It was a newly emerging region 
which developed to a large AR having f3^8 mag- 
netic complexity during its rapid evolution. It con- 
sisted of four large regions of opposite polarities 
in quadrupolar configuration. Figure QJleft col- 
umn) shows the evolution of NOAA 11158 during 
2011 February 13-16 in HMI intensity maps. The 
prominent positive polarity sunspots are labeled 
as SP1,SP2,SP3 and the negative polarity spots 
as SN1,SN2, SN3 for identification. LOS contours 
are overlaid on the intensity image showing the 
respective polarity distribution. 

The spatial evolution of the AR shows a large 
shearing motion of SP2 that rotated around SN2 
about its umbral axis during 2011 February 13-15. 
It then detached and moved towards SP3 along 
with small patches of both polarities appearing 
and disappearing over short periods of time. This 
motion appeared to have created a twist in mag- 
netic fieldlines connecting these spots. A careful 
examination of the animation made from mag- 
netograms and intensity maps revealed a signifi- 
cant counter clock-wise rotation of SN1 during the 
same period, while a small positive polarity region 
SP1 located to the north of SN1 rotated in the 
counter clock-wise direction along with a proper 
motion away from SN1. The rotations of SN1 and 
SP1 increased the twist of the field lines, and the 
magnetic non-p o tentia lity of the sigmoid structure 
(|Canfield et all (|l999l) l ' Several mass expulsions 
were launched intermittently from this region, as 
seen from the quick look images in AIA. These 



turned into CMEs as confirmed by STEREO ob- 
servations. 

In order to inspect the magnetic complexity or 
twist contributed by the observed shearing mo- 
tions of the magnetic foot points, we computed 
the helicity injection rates using the temporal se- 
quence of magnetograms of the AR. Figure Upright 
panel) shows the computed helicity density maps 
corresponding to the HMI continuum intensity im- 
ages (left panel). The dark (white) patches in the 
right panel represent negative (positive) helicity 
density according to the usual convention. Con- 
tours of LOS magnetic field at [-150, 150]G levels 
are overlaid for a better visualization of helicity 
density with respect to the magnetic polarity. Ev- 
idently, negative polarity region of SN1 injected 
negative (dark) helicity during 2011 February 14- 
15 which is also consistent with its physical CCW 
rotation. On the other hand, SP2, SN2 and SN3 
injected positive (white) helicity along with neg- 
ative (dark) helicity in some small patches. We 
expect that the nature of motions in these areas 
could have influenced the helicity pattern there. 
One can quantitatively estimate these motions by 
means of the twist or structural complexity of 
these regions. 

The photospheric maps of helicity (and its in- 
jection rate) provides spatial information about 
the basic properties of a link between the activ- 
ity and its sub-photospheri c roots as reflected by 
the flux emergence process. Ueong fc Chad <l2007h 
have shown that helicity is mostly injected while 
fluxes emerged in the AR suggesting it to be the 
major source of helicity injection. Flux cancela- 
tion process, on the other hand, resulted in a loss 
of coronal magnetic helicity, or inverse helicity in- 
jection. We thus infer that the AR possessed two 
main sites, Rl and R2, of unstable energy stor- 
age systems marked by the rectangular boxes in 
Figure [TJ These sites had distinctly different in- 
jection of helicity density corresponding to the flux 
(or foot-point) motions, polarities and activity. 

Figure [2] shows transverse velocities in the rect- 
angular regions Rl (top row) and R2 (bottom 
row) of NOAA 11158 overlaid on the correspond- 
ing maps of helicity density during the X2.2 flare 
on February 15/01:44UT. Also overlaid are the 
contours of the LOS magnetic field. Maximum 
rms velocities in the range of 0.6-0.7 km s _1 were 
found in the AR. This map indicates the highly 
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Fig. 1. — (Left column) The daily HMI continuum intensity maps of AR NOAA 11158, and (Right column) 
the corresponding helicity density maps (scaled to ±0.05 x 10 20 Mx 2 cm~ 2 s _1 ) computed from Eq|H The 
field of view is 275 x 200 arcsec 2 . The overlaid gray and black contours correspond to LOS magnetic fields at 
[-150,150]G levels, respectively. Rectangular boxes in intensity image of 2011 February 14 mark the selected 
sub-areas Rl and R2 in which velocity flows are shown in the subsequent figures. 
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Fig. 2. — Transverse velocity field vectors as inferred from DAVE technique superposed on helicity density 
maps with the LOS magnetic field contours for the rectangular regions of Figure [I] - Rl (Left column) and 
R2 (Right column). Spiral velocity patterns in (b) indicate vortical like motion in the sunspot penumbra 
due to rotation. Sites of negative helicity injection are seen along magnetic polarity inversion line in (e)-(h), 
particularly as marked by a rectangular box in (g). SP1 site is encircled in (b). 
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dynamic state of the AR by a factor of two higher 



than that reported earlier (jKusano et al.l [2002) 



This may have to do with the higher sensitivity 
and spatial resolution of the HMI instrument. 

Spiral or vortex like velocity patterns are obvi- 
ously related to the counter rotation of SN1 and 
the CCW rotating, emerging region SP1 of en- 
circled in Figure EJb) . A noticeable evidence is 
that the sub-region Rl possessed negative helicity 
density distribution which is consistent with the 
chirality associated with the physically observed 
counter rotation of SN1. Because of the continued 
shearing motions at the interface of SP2 and SN2, 
the flow field vectors almost aligned with the po- 
larity inversion line (PIL) as seen in (e-g). Inter- 
action of fluxes with this shear motion can squeeze 
and converge the flux in both SP2 and SN2. The 
field lines were stressed and twisted by this mo- 
tion leading to the storage of free energy adequate 
to account for the release in the energetic X2.2 
flare of February 15/01:44UT. By examining the 
spatial distribution of helicity change in nearly all 
the flare events, we noticed sudden appearance of 
negative patch of helicity in the regions of posi- 
tive helicity. This process is particularly well dis- 
cernible in panels (g)-(h) in the region marked by 
the rectangular box that was the site of the ener- 
getic X2.2 flare. 

We have plotted the disk integrated GOES soft 
X-ray flux (1-8 A channel) during February 11-17 
in Figure E^top) where the peak times of flares 
of NOAA 11158 are marked by arrows. After 
its birth, the AR gradually evolved during 2011 
February 11 - 13 as evident from the monotonic 
increase of fluxes in both polarities corresponding 
to 3 x 10 21 Mx (Figured middle). Then followed 
a rapid phase of flux emergence (of 9 x 10 21 Mx) 
during February 13-14 after which it reached a 
plateau. Also plotted is the flux imbalance, i.e, 
the ratio of the net flux and absolute total flux 
in the AR. The dominance of negative flux dur- 
ing February 13-15, and thereafter of the positive 
flux is evident. Flux variations occurred in the 
range of (9.5-12.5) x 10 21 Mx with the unbalance 
within 10%. A significant flux decrease occurred 
in both polarities by ~ 1 x 10 21 Mx at the time 
of the X2.2 flare. This implies that flux amount- 
ing to ~ 1 x 10 21 Mx might have been annihilated 
during this energetic flare. The unusual rotating 
sunspots along with the increased fluxes indicated 



emergence of highly twisted fluxes fr om the sub- 
photospheric region (jLeka et al ] ll99fih . and not re- 
sulting from the surface flows alone. Most of the 
flare and CME activity of this AR occurred only 
after February 13/12:00UT, indicating that the 
rapid flux emergence could have played important 
role in triggering the transients. 

In Figure [^bottom), we have plotted the time 
profile of helicity injection rate, which is the sum- 
mation of helicity density over the AR divided by 
the time interval. Also plotted is the accumu- 
lated helicity, i.e, the integrated helicity change 
rate over time. The total accumulated helicity is 
estimated as 13.30 xl0 42 Mx 2 during the six day 
period of 2011 February 11-17, with the peak he- 
licity rate of 32xl0 40 Mx 2 h~ 1 . The occurrence 
times of the CMEs associated with the AR are 
marked by arrows in this panel for reference. An 
impulsive variation of helicity injection rate due 
to injection of negative helicity is discernible dur- 
ing the X2.2 flare and the concomitant CME. The 
helicity injection rate decreased during the period 
February 14/ll:00-February 15/13:00 UT, and in- 
creased thereafter till February 16 along with fluc- 
tuations in the range 2-4 x 10 40 Mx 2 h _1 . 

The shaded vertical bars in Figure [3] corre- 
spond to the impulsive variations at the time of 
the flares. Some of the flares were also associated 
with CMEs or plasma ejections. Due to the fluc- 
tuations, it is not possible to attribute the varia- 
tions associated with small flares, therefore they 
are not marked by shaded bars. Further, we ex- 
amined the spatial distribution of helicity change 
rate around the time of the X2.2 flare as plotted 
in Figure [21 Due to the appearance of negative 
sign patch of helicity change rate in the region 
of dominant helicity of opposite sign(panel f), an 
abrupt decrease of the positive injection rate oc- 
curred on February 15/02:00UT that was appar- 
ently seen in the time evolution plot. Similar fea- 
tures were also found for the M6.6 flare on Feb 
13/17:28UT, C7.7 flare on Feb 14/13:47UT and 
M1.6 flare on Feb 16/14:19UT events(see Figure 
10-12). However, there are also other such impul- 
sive variations which may be related to variations 
at sites other than the flare sites. This may arise 
as the profiles are obtained by integrating over the 
entire AR. It is important to examine these varia- 
tions by theoretical arguments to explain any pos- 
sible role they might have played in the initiation 
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AR 11158 




Fig. 3. — Top: Solar disk integrated GOES Soft X-ray flux during February 11-17, 2011. The arrows on top 
panel indicate the start times of flares in AR NOAA 11158. Middle: Time profiles of the magnetic fluxes 
and flux imbalance in the AR. Bottom: The computed helicity rates integrated over the whole AR. Arrows 
in this panel indicate the onset time of CMEs that were launched from this AR. Shaded vertical bars show 
the flares (and associated CMEs) with sharp variation of helicity injection rates. 
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Fig. 4. — (Left column) The daily HMI continuum intensity maps of AR NOAA 11166, and (Right column) 
the corresponding helicity density maps computed from Eq The field of view is 350 x 200 arcsec 2 . 
The overlaid gray and black contours correspond to LOS magnetic fields at [-150,150]G levels, respectively. 
Rectangular boxes in intensity image of March 9 mjirk the selected sub-areas in which velocity flows are 
shown in the next figure. Emerging fluxes from sunspot periphery are denoted as FN2 and FP3 on March 
11/22:00UT. 



of eruptive events. 

3.2. AR NOAA 11166 

AR NOAA 11166 appeared on the east limb 
of solar disk on 2011 March 03 at the location 
N10E64. We monitored its activity during the pe- 
riod of 2011 March 6-11 in which it produced a 
large XI. 5 flare, two M-class flares and several C- 
class flares, some of which were also associated 
with plasmoid ejections or CMEs. Table Q] lists 
the flares and CMEs of this AR. Daily evolution 
of the AR in the period of March 8-11, 2011 is 
shown in Figure rjjleft column). 

The major sunspots of the AR are labeled as 
SP1, SP2, SN1 and SN2. The identification of SP2 
was somewhat unclear before March 10 as several 
small umbrae were spread over its location. They 
moved and coalesced to form SP2. Polarities of the 
respective sunspots are identified by the overlaid 
LOS magnetic field contours. This AR also con- 
sisted of a complex magnetic configuration with 
two positive (SP1, SP2) and two negative (SN1, 
SN2) polarity sunspots located within the sur- 
rounding diffused fluxes. Emerging and moving 
flux regions, FP3 and FN2, were identified in the 
course of the evolution in the sunspot periphery 
(March 11/22:00UT panel), having opposite sign 
to that of their native sunspots. However, there 
were no intrinsic rotating sunspots or flux patches 
as observed in the case of AR NOAA 11158. 

We computed helicity density for AR NOAA 
11166 during its evolution in the period 2011 
March 6—11. The corresponding maps are shown 
in Figure fright column) . Locations of helicity 
density of mixed sign were distributed all over 
the AR. The peripheral sites of the sunspots ex- 
hibited helicity density predominantly of negative 
sign. However, patches of negative sign helicity 
were also observed embedded in the positive helic- 
ity site of the flare (March 09/23:00UT panel). For 
further examination, we consider two sub-areas Rl 
and R2, as marked by the boxes in this panel. 

Horizontal, or transverse, velocity vectors cor- 
responding to the tracked flux motions are plot- 
ted in Figure [5] separately for Rl (top row) and 
R2(bottom row) during the time of the XI. 5 flare. 
The rms velocities are found to have the maxi- 
mum values in the range 0.5-0.7 km s . Strong 
moat flows were systematically dominant in both 



regions from the peripheral regions of sunspots in 
addition to the shear flows. Persistent strong shear 
motions due to the merging SP2 group were iden- 
tified in R2. These flows appear to collide head 
on with those from SP1 resulting in the submer- 
gence. Flux emerging region was identified from 
the diverging flow field shown in the encircled re- 
gion in (c). From this region, flux moved towards 
R2 as the AR evolved. Negatively signed helic- 
ity is prominent in (d)-(f) of R2 as marked in the 
rectangular box in (e). Negative helicity injected 
mostly from SP2 at the emerging flux site FN2 
enclosed in the rectangular boxes in (a)-(d) in the 
sub-region Rl. The intense negative helicity injec- 
tion is evident from (e)-(f) at the time of the large 
CME that was followed by XI. 5 flare on March 
09/23T3UT. 

The disk integrated GOES soft X-ray flux (1- 
8A channel) during 2011 March 6-12 is plotted 
in Figure [Hltop). The arrows in this panel in- 
dicate the peak time of flares in NOAA 11166. 
During the disk transit of the AR, fluxes of both 
polarities increased corresponding to 5 x 10 21 Mx, 
with the imbalance varying below 6% (Figure [6j 
middle). As observed for NOAA 11158, a rapid 
flux emergence occurred in this AR too during 
March 7-9. Thereafter, only small variations asso- 
ciated with local cancelations/emergence of about 
~ 1 x 10 21 Mx took place. Positive flux dominated 
in the AR during March 7-11, and then a near 
balance was established. 

It is worth noticing that magnetic fluxes in both 
polarities decreased by ~ 0.9 x 10 21 Mx at the time 
of the XI. 5 flare. The net flux imbalance, increas- 
ing prior to the flare, reduced significantly after 
the flare. This may indicate that the flux amount- 
ing to ~ 1.8 x 10 21 Mx might have been annihilated 
during this flare. Most of the flares and CME ac- 
tivity of this AR occurred only after March 8, sug- 
gesting that the rapid emergence of fluxes could 
be an important factor for triggering these tran- 
sients. This result is similar to the case of AR 
NOAA 11158. 

Temporal evolution of helicity injection rate 
and the accumulated helicity for NOAA 11166 are 
shown in Figure [BJbottom) with arrows marking 
the times of the CMEs. As expected, these pa- 
rameters increased in the first phase correspond- 
ing to the flux eme rgence, in agreement with 
Jeong fc Chad (|2007l ) that helicity is mostly in- 
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Fig. 5. — Transverse velocity field vectors in the rectangular region Rl (Left column) and R2 (Right column) 
of Figure|4]as inferred from DAVE overlaid on the helicity density maps(scaled at ±0.05 x 10 20 Mx 2 cm~ 2 s~ 1 ) 
with iso-contours of LOS magnetic field during the XI. 5 flare. The encircled region shows diverging flow. 
Regions of negative helicity are seen in the rectangular boxes. 
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6.— Same as Figure[3]but for AR NOAA 11166. 
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jected while the fluxes emerged. Total helicity ac- 
cumulated during the six days' period of the AR's 
evolution was estimated to « 9.5 x 10 42 Mx 2 . The 
maximum helicity injection occurred during 2011 
March 8 at the rate of (30 x lO^Mx 2 ^ 1 ). There- 
after, it reduced gradually to the minimum rate at 
-10 x lO^Mx 2 /^ 1 as observed on 2011 March 10. 
The coronal helicity of the AR is likely to be pos- 
itive as a result of the positive helicity injection. 

Vertical shaded bars in the figure indicate the 
epochs of sudden variation of helicity rate at the 
onset of the major flares of NOAA 11166. Some 
of these flares were also associated with CMEs. 
For instance, on 2011 March 9/21:45UT, a fast 
CME was followed by the XI. 5 flare during which 
a large dip in helicity rate occurred. The same 
is the case on March 10 with M-class flare. This 
could be an important finding of this study that 
the injection of opposite helicity led to initiation 
of CMEs followed by flares. 

3.3. Dependence of Helicity Injection Rate 
on the DAVE Parameters 

Computation of helicity injection rate involves 
the measurement of magnetic field and the in- 
ferred horizontal velocities. Apart from the errors 
in the measurements, the computations involving 
the DAVE method for deriving velocities depend 
on two main parameters viz., the time interval be- 
tween two successive magnetic maps, At, and the 
DAVE window size. For obtaining optimized re- 
sults, horizontal displacements of features during 
the time interval At should be large enough to be 
well determined by DAVE. Also, these displace- 
ments should be smaller than the selected win- 
dow size. To check our results for consistency, we 
carried out the DAVE calculations using the time 
intervals At = 12, 24 and 36 minutes, while, keep- 
ing the window size fixed at 21 x 18 pixels. Then, 
calculations were carried out for different window 
sizes, viz., (21x18, 15x12, 9x6) while keeping At 
fixed at 12 minutes. Furthermore, to avoid the 
effect arising from noise, we used a threshold of 
magnetic field at 10G, which is the HMI preci- 
sion. As the HMI provides 12 minute averaged 
data products, we averaged them corresponding to 
our calculations at 24 (2 maps) and 36 (3 maps) 
minutes. 

The dependence of helicity injection rates on 
time interval At is shown in Figure [TJtop row) for 



NOAA 11158. The scattered data are fitted by 
straight line in the least square sense. Due to the 
large volume of data, this computation is tedious 
and time consuming. Therefore, results are shown 
here only for NOAA 11158, however, they are also 
valid for other ARs observed by the HMI. There is 
an additional issue of unequally spaced data points 
to be addressed in case, for example, we intend to 
plot the results for At=36 with At = 24 minutes. 
For such cases, we used a cubic spline interpolation 
(cf. lPress et al. (1992)), to get corresponding ab- 
scissa values for the ordinate points or vice- versa. 
Essentially, this algorithm employs cubic polyno- 
mial between each pair of data points with the 
constraint that the second and first derivative of 
that polynomial is same at the end points so that 
resulting values are smooth. 

Table 2: Helicity injection rates at different DAVE 
parameters 



DAVE parameters 


AR 11158 


At 


Window size 


dH/dt 


Acc Hel 


min 


pix 2 


min 


max 




12 


21x18 


-19.13 


32.47 


13.30 


24 


21x18 


-8.14 


26.90 


12.52 


36 


21x18 


-1.06 


22.51 


12.15 


12 


21x18 


-19.13 


32.47 


13.30 


12 


15x12 


-18.32 


33.76 


14.83 


12 


9x6 


-15.81 


37.09 


16.14 



Units of dH/dt are 10 40 Mx 2 h- 1 and Accu Hel 
are 10 42 Mx 2 

It can be observed from the scatter plots that 
the helicity rates decrease slightly as the time in- 
terval At is increased from 12 min to 36 min. This 
implies that short-lived features contribute to he- 
licity rates but their effects are rather small. In 
addition, best slope (0.91) as well as correlation 
(0.96) exists between the time intervals of 36 and 
24 minutes. The variation in accumulated helic- 
ity is less than 9%. This implies that averaging 
in time has negligible effect on the results and 
the DAVE method is able to reproduce nearly the 
same results. 

The dependence of helicity injection rate on 
window size by keeping the time interval At fixed 
at 12 minutes is shown in Figure [7l A significant 
dependence exists on the window size with rea- 
sonable correlation coefficients but a wide range 
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Fig. 7. — Dependence of helicity injection rate for AR NOAA 11158 on (Top row) the time interval At, and 
(bottom row) the window size. 



1G 




Fig. 8. — (Top) Helicity injection rate, and (bottom) the accumulated helicity obtained for the rectangular 
regions Rl and R2 of Figure [Qfor AR 11158 (Left) and FigureHfor AR 11168 (Right). 



of slopes. Accumulated helicity increased for the 
DAVE window size from 21x18 to 9x6. A larger 
scatter or lower correlation coefficient (Corr = 
0.97) exists for window size 21 x 18 with 9x6 pixel. 
This is expected, because affine velocity profile 
is not valid for large window size. Also, DAVE 
method estimates the velocity fields as constant 
within the window as does the LCT method. A 
total variation of 21% is found, however, with the 
same trend of helicity injection rate profiles which 
is discernible in correlation coefficient in the plot. 
This dependence is o pposite of the earlier report 
by iTian et ah ( 2011 ) . A maximum velocity of 1 
km-s -1 during the time interval of 12 min cor- 
responds to a plasma displacement of an arc-sec. 
Hence, for the window size of 4".5x3"(9x6 pix), 
the issue of features overflowing out of the window 
should not pose problem. 

Table [2] lists the minimum and maximum val- 
ues of helicity injection rates (dH/dt, in units of 
10 40 Mx 2 h _1 ) and the accumulated helicity (Acc 
Hel, in units of 10 42 Mx 2 ) for the computational 
runs carried out with various DAVE parameters 
as mentioned above. 

We, thus, found the measured helicity injec- 
tion rate to depend on the time interval between 
the two successive magnetograms, i.e., the obser- 
vational cadence. The selected window size, i.e., 



the spatial resolution, also influences the measured 
quantities. Our analysis suggests that it is better 
to use images averaged over up to 24 minutes with 
relatively small DAVE window size subjected to 
the overflow condition as mentioned above. These 
are important considerations to derive reasonable 
and meaningful results in addition to optimizing 
the computations involving large data-sets. 

4. Discussions 

Free energy storage and release are some of the 
most important problems in the eruption physics 
of the Sun. There are essentially two effects that 
can supply magnetic free energy and helicity from 
below the solar surface to the corona. Flux emer- 
gence is the process in which vertical motions carry 
magnetic fluxes through the photosphere. If the 
sub-surface fluxes emerging through the photo- 
sphere are already twisted, then it will contribute 
to the injection of helicity (c.f., the 1st term in 
Eq [J). Computation of this term requires the 
knowledge of the vertical component of velocity 
and the horizontal or transverse component of 
magnetic field. Flux motions in the form of ro- 
tation or proper motions are another process that 
may efficiently supply helicity injection (c.f., the 
2nd term in Eq[5]). Solar differential rotation con- 
tribution to helicity injection is rather small, less 
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Fig. 9. — Accumulated helicity with total absolute flux computed for NOAA 11158(Left) and NOAA 
11166(Right) . Events were labeled and shown by circles on the plot. 



than 10% of that contributed by the flu x motions 
(|Chae et alj|2004 iDemoulin et "al.ll2002h . and has 
only a much long te rm effect on helicity accumu- 
lation (jDeVorell2000h . 

Magnetic helicity is a physical quantity having 
a positive or negative sign, representing a right- 
handed or left-handed linkage of magnetic fluxes, 
respectively. This means that if positive and neg- 
ative helicities co-exist in a single domain, mag- 
netic reconnection can cancel magnetic helicity by 
merging magnetic flux systems of opposite helic- 
ities. Our computations of magnetic helicity in 
both ARs revealed that the distribution of helicity 
injection is highly complicated in time and space. 
Even the sign of helicity injection often changed 
within the AR. For instance, the sub-area Rl in 
AR 11158 of negative sign helicity produced most 
of the CMEs. On the other hand, most of the 
flares originated in R2 of positive helicity, includ- 
ing the X2.2 flare. In AR 11166, however, helicities 
of mixed sign were observed in both sub-areas Rl 
and R2, which produced flares and CMEs. Specif- 
ically, spatial distribution of helicity injection at 
the time of XI. 5 flare is shown in Figured] 

We have plotted the total injected quantities 
for R1-R2 of the respective ARs in Figure [5] In- 
jection of helicity in a region of dominant opposite 



sign can be understood as a sudden dip in the time 
profile plot. Of course, the corresponding spatial 
information is lost in the averaged quantity. The 
advantage of using localized analysis of selected 
sub-areas in the ARs is that it reduces complex 
variations occurring over a much larger area while 
showing only the variations occurring in the areas- 
of-intcrcst. It also enhances the dips correspond- 
ing to the identified events marked by the arrows. 
However, it is important to identify the location 
of individual event in order to correctly attribute 
a particular change of helicity rate to it. NOAA 
11158 was essentially a positive helicity injecting 
region, while its sub-region Rl had a negative in- 
jection rate and accumulated quantity due to pres- 
ence of rotational motion. We expect that as the 
sunspots SN1 and SP1 rotated, the injection rate 
increased to a maximum rate of -16 x lO^MA -1 
on February 14/18:00UT. A total helicity accu- 
mulation of —4.3 x 10 43 Mx 2 occurred during the 
six day period from this region. Noticeably, a 
steep accumulation was occurred during Feb 14- 
15 along with many observed expulsions shown by 
arrows. This may be interpreted as shedding of 
excess helicity from the corona in the form of erup- 
tive events. The steep accumulation of helicity by 
monotonic injection rate, therefore, is suggested 
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to a signature of expulsions. Accumulated helic- 
ity amounts to 11.5 x 10 43 Mx 2 from sub-region R2 
and from Feb 13 onwards, steep accumulation ob- 
served which could be mostly associated with ob- 
served large shear motion of SP2. In NOAA 11166, 
Rl had positive injection rate with positive accu- 
mulated helicity of 11 x 10 40 Mx 2 h _1 . However, a 
varying injection rates of both signs with vary ac- 
cumulated helicity were evident in sub-region R2. 
This trend can not be interpreted correctly unless 
flux motions entirely localized in the considered 
sub-region. A prominent peak of helicity rate was 
found around the time of the XI. 5 flare in R2 of 
AR 11166 that was not obvious in Figure middle 
panel). 

It has been suggested earlier by several work- 
ers that magnetic helicity must play an impor- 
tant role in flares as a substantial amount of he- 
licity accumulation is found before many events 
dKusano et al.l(ll995l ); [Kusano k Nishikawal (| 19961) ; 
iKusano et al.l (|2002f) ). However, the correlation 
between various magnetic field parameters and 
the flare index of an AR is not high irrespective of 
the method used. This is an intrinsic problem for 
flare forecasting as the occurrence of a flare de- 
pends not only on the amount of magnetic energy 
stored in an AR, but also on how it is triggered. 
Thus, it appears that helicity accumulation might 
be a necessary, but insufficient condition for the 
flares requiring a trigger even if a magnetic sys- 
tem has enough non -potentiality. For instance, 
Kus ano et al. (2002) suggested that coexistence 



of positive and negative helicities is crucial for the 
onset of flares. 

Careful three-dim ensional simulation s have 
been carried out by iLinton et al. ( 2001 ) to ex- 
plore the physics of flux tube interaction for the 
co-helicity (same sign) or counter-helicity (oppo- 
site sign). According to them, counter-helicity 
presented the most energetic type of slingshot in- 
teraction in which flux is annihilated and twist is 
canceled. On the other hand, co-helicity exhibited 
very little interaction, and the flux tubes bounced 
off resulting in negligible magnetic energy release. 

Magnetic helicity in the solar corona is closely 
related to the photospheric magnetic shear, which 
is usually defined as the extent of alignment of 
the transverse component of magnetic field along 
the neutral or polarity inversion line (PIL). Based 
on this idea, IKusano et al.l (|2004l ) performed a nu- 



merical simulation by applying a slow foot point 
motion. This motion can reverse the preloaded 
magnetic shear at the PIL resulting in a large scale 
eruption of the magnetic arcade through a series 
of two different kinds of magnetic reconnections. 
They proposed a model for solar flares in which 
magnetic reconnection converts oppositely sheared 
field into shear free fields. 

We have found evidence of such interactions 
during the X-class flares and associated CMEs as 
seen in Figure [2] in the form of continued shear- 
ing motion of SP2 around SN2 in AR 11158. 
Similar motions are also associated with SP2 in 
AR 11166. Correspondingly, sudden impulsive 
peaks appeared in the profiles of helicity injec- 
tion(shaded bars in Figure |3l [HJ due to the in- 
jection of negative signed helicity. Opposite he- 
licity flux tubes can interact easily leading to re- 
connection, thereby unleashing explosive release of 
magnetic energy. Impulsive variations of magnetic 
helicity injection rate associated with eruptive X- 
and M- class flare s accompanied with CMEs were 
reported a l so by Moon et al.l (|2002j). Recently, 
iPark et al.l (|2010al ) conjecture that the occurrence 
of the X3.4 flare on 2006 December 13 was in- 
volved with the positive helicity injection into an 
existing system of negative helicity. Further, so- 
lar eruption triggered by opposite helicity o f two 
magnetic flux systems ([Romano et al.l 2011 ) , and 
occurrence of flares in rel ation to spatial distri- 
butio n of helicity density (|Romano fc Zuccarello 



20111 ) were reported in the literature. The main 



drawback of these findings is that the time span 
between two magnetograms is more than duration 
of the flare (> 96m), so time rate of helicity could 
not be easily resolved at the onset time of the flare. 
Thus, our results further support the idea that op- 
posite helicity flux tubes reconnect to annihilate 
magnetic fields. 

Using the above mentioned theoretical argu- 
ments, we infer from our observations that the 
sudden impulsive variations in helicity injection 
are due to the mechanism by which the highly 
sheared field configuration caused by photospheric 
motions is converted to shear-free field by releas- 
ing excess energy. In both the ARs we examined, 
minimum injection rates are found after X-flarc 
associated with CME (after/before flare), indicat- 
ing the condition of lower shear. Fluctuations in 
the temporal profiles of helicity injection rate at 
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times other than shown in the figure include sys- 
tematic as well as computational errors. These 
errors dilute the exact change of the signal, there- 
fore small magnitude flares did not show identifi- 
able variations. Area-averaging over small regions 
is expected to reduce such fluctuations as in Fig- 
ure HI However, it is essential to correctly identify 
location of individual event within the AR. Fur- 
ther, our analysis needs to be suitably rechecked 
by using vector magnetic fields and the emergence 
term. Once this is confirmed, one may be able 
to use the helicity injection rate as the forecast 
system for eruptive events. 

In addition, there are concerns about the nature 
of changes observed in magnetic field measure- 
ments during major energetic transients. Tran- 
sient magnetic field and Doppler velocity changes 
could arise due to spectral line profile changes from 
absorption to emission thermal and non-thermal 
process affecti ng the photospheri c layer during 
X- cla ss flares ( Qiu fc Gary 20031 : Maurva et al 
20111 ). However, these transient effects are di- 



luted if time averages of 12 min or longer are used 
since transient variations typically occur during 
the short period of the impulsive phase of the flare, 
i.e., about 3-4 minutes. In order to overcome this 
difficulty, one may choose to ignore any impulsive 
helicity variations if At is within this limit. 

Our computed helicity rates involving photo- 
spheric flux motions include the flux em e rgenc e 
term as explained by iDemoulin fc Berger (2003). 
By a simple geometrical argument, horizontal 
foot-point velocity (u, here DAVE velocity) can be 
written in terms of horizontal and vertical plasma 
velocities, v^, v n , respectively: 



Vh 



-Bv 



(5) 



From this relation, it is not possible to infer as 
to which term, flux emergence or flux motions, 
governs the level of activity of the ARs. To resolve 
this difficulty, we plot the integrated absolute flux 
and accumulated helicity computed over the AR, 
as shown in Figure El 

Evidently, accumulated helicity increased mono- 
tonically with the emergence of magnetic flux in 
the AR in its emerging phase (marked by the verti- 
cal dashed line for NOAA 11158). After this phase 
followed the next, the active phase, where an ap- 
preciable increase of helicity occurred with only 



small variation in the flux. This rapid increase in 
helicity could be interpreted as the dominant con- 
tribution of the flux motions. Intermittent mass 
expulsions in the form of CMEs transferred away 
the excess helicity. The extent of this transfer, 
however, is not clear from this plot, although one 
can make plausible conclusions from the timings 
of the flares and CMEs. 

It is evident that flux emergence is necessary 
and their motions are crucial for increasing the 
complexity of magnetic structures in the form of 
shear and helicity. Therefore, these motions addi- 
tionally accounted for the accumulated amount of 
emergence term. However, as mentioned earlier, 
more careful and exact calculations are required 
to implement the emergence term by using vector 
magnetic field data. We intend to carry out such 
improved calculations in the future when the re- 
quired vector magnetograms become available in 
order to further confirm our results. 

5. Summary 

We have studied the evolution of magnetic 
fluxes, horizontal flux motions, helicity injection 
and their relationship with the eruptive transient 
events in two recent flare (CME) productive ARs, 
NOAA 11158 and NOAA 11166 of 2011 Febru- 
ary and March, respectively. We have used high 
resolution, high cadence data provided by SDO- 
HMI for these ARs which were in their emerging 
and active phases. The emerging ARs consisted 
of rotating sunspots with increasing flux indicating 
twisted flux emergence from the sub-photospheric, 
convection zone. This indicated the transfer of 
twist or helicity injection through the photosphere 
to the outer atmosphere. The results are summa- 
rized as follows: 

1. Strong shear motions that include rotational 
and proper motions played significant role 
in most of the events in addition to the 
flux emergence. Such motions are crucial in 
twisting or shearing the magnetic field lines 
and for further flux interactions. 

2. AR NOAA 11158 consisted of a CME-prone 
site of rotating main sunspot along with 
emerging flux of opposite sign and counter 
rotating moving magnetic feature. It also 
had a flare-prone site consisting of self- 
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rotating sunspot moving about a sunspot 
of opposite sign, leading to flux interaction. 
Both these sites led to the formation of sig- 
moidal structures, which are unstable, and 
more likely to produce eruptive events. A 
strong CME on 2011 February 14/17:35UT 
occurred in the former site and a white 
light, energetic X2.2 flare on 2011 Febru- 
ary 15/01 :44UT occurred in the later site. 

3. AR NOAA 11166 was in its active phase 
with further increasing content of flux. 
Group motions of diffused fluxes merging 
to form bigger sunspot manifested major 
shear motions in addition to outward flows 
from sunspot. A large CME on 2011 March 
09/21:45 followed by an XI. 5 flare was one 
of the major events in this AR as a result of 
these motions. 

4. We found spatial and temporal correspon- 
dence between flaring sites and negative 
hclicity injecting regions. Area integrated 
temporal profiles showed sudden variations 
in helicity injection rates at the onset of 
flares or CMEs. This supports the valid- 
ity of triggeri n g mec hanism as proposed by 



Kusanoetal 



(120041 ) . In both the studied 
ARs, we found minimum injection rates af- 
ter X-flare with CME (after/before flare), 
indicating shear free conditions. A further 
improvement of these investigations requires 
vector magnetic field data in order to derive 
the LOS velocity for computing the emer- 
gence term of helicity injection. 

AR NOAA 11158 injected 13.30 x 10 42 Mx 2 
helicity while AR NOAA 11166 injected 
9.5 x 10 42 Mx 2 helicity during the six days' 
period of their evolution. These are consis- 
tent with the previou sly reported range of 
helicity injection (eg., Park et al. ( 2010bl )). 
It appears that due to the presence of rota- 
tional motion, the former AR accumulated 
larger amount of helicity accounting for its 
greater activity in the form of flares and 
CMEs. 



6. Apart from the instrumental and computa- 
tional errors, the estimation of helicity in- 
jection rates also depends on DAVE param- 
eters used to track the motion of the fluxes. 



This is evident from a 9% decreases in the 
accumulated helicity by increasing the time 
interval between successive magnetograms, 
or the cadence. Similarly, it increased by 
21% when smaller window sizes were used 
in the computations keeping the time inter- 
val fixed. 

We plan to improve the results obtained in this 
study by including the emergence term and fur- 
ther establish the relation of helicity injection with 
eruptive events. This is expected to be carried out 
with the availability of the standard HMI vector 
magnetic field maps. It is also required to study a 
larger sample of ARs to generalize our understand- 
ing of the role of helicity injection in the onset of 
the eruptive events. 

The data have been used here courtesy of 
NASA/SDO and HMI science team. We thank 
Dr. Etienne Pariat for checking our helicity pro- 
gram with comments and suggestions. Special 
thanks to Mr.Jigar Raval, Mr.Anish Parwage for 
their help in running program on one of the nodes 
of HPC cluster at PRL computer center. 
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Fig. 10.— M6.6 flare on Feb 13 at 17:28UT. Top row: Rl and Bottom row: R2 




Fig. 11 — C7.7 flare on Feb 14 at 13:47UT. Top row: Rl and Bottom row: R2 
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Fig. 12.— M1.6 flare on Feb 16 at 14:19UT. Top row: Rl and Bottom row: R2 
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